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Abstract 

For a general spin glass model with asymmetric couplings we prove a family of 
identities involving expectations of generalized overlaps and magnetizations in the 
quenched state. Those identities holds pointwise in the Nishimori line and are 
reached at the rate of the inverse volume while, in the general case, they can be 
proved in integral average. 
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1 Introduction and results 



Overlap identities have played and continue to play a central role in spin glass statistical 
mechanics since the appearance of the Parisi solution [MPV] of the Sherrington Kirk- 
patrick model [SK] . The replica symmetry breaking theory contains indeed, as a built-in 
ansatz property, a family of identity for the overlap expectations that have been since 
then classified as replica equivalence and ultrametricity factorization properties. While 
the first is now largely understood for both mean field and short range finite dimensional 
models [CGi, CGi2] the second is still an open conjecture even within the mean field 
cases. The ideas to obtain the rigorous proof of the identities trace back to the papers 
[AC, GG] where the invariance property of stochastic stabihty and the role of the en- 
ergy fluctuations for the spin glass quenched measure were introduced (see also [Ba] for 
a different and original derivation). The mentioned properties state a striking feature for 
the quenched measure on overlap expectations: the overlap moments do obey sum rules 
that reduce the distributional degrees of freedom and, at least in the mean field case, it 
is expected that the entire distribution be identified by the simple overlap distribution. 
Those relations are expected to hold everywhere but on isolated singularities and they 
can be in fact rigorously proved in f3 Riemann integral average [CGi2, T, B]. 

We present in this paper a new family of identities in terms of generalized overlap 
and magnetization that hold when the interactions are not centered which are proved in 
integral average. Moreover we find the remarkable result that in the Nishimori line they 
hold everywhere with respect to the parameters and are reached at the expected rate of 
the inverse volume. 

We can illustrate the results proved in this paper considering the Edwards- Anderson 
spin glass model with non symmetric Gaussian interactions of variance and non zero 
average /i. Considering the quenched measure over multiple copies of the system subject 
to the same disorder (for a precise definition see Section 2) the link overlaps and the 
link magnetizations fulfill an infinite family of identities when integrated in the inverse 
temperature /3 over arbitrary intervals. Those identities at the lowest moment read, for 
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all /X and in the thermodynamic limit: 



I dp \ < -2gi 2 + 891,292,3 - 6^1,293,4 > 
+ P < 4gi,2m3 - 4gi,2"^i > 

+ < ml — mim2 > | = (1-1) 



/•ft 

/ dp I P^{< > -6 < gi,2 >< 92,3 > +6 < gi,293,4 > - < 9l,2 >^) 

+ P{2 < gi,2mi > -4 < 5i,2m3 > +2 < 51,2 >< mi >) 

+ < mim2 > - < mi >^ | = (1.2) 

Since we prove that the magnetization has vanishing fluctuation in //-average in the infinite 
volume limit 

/ dii{< m^ > - < m >^) = (1.3) 
J Ml 

then, when integrated over arbitrary boxes in the (/3, //) plane, the previous identities 

become 

dfx dp\P^ < -2ql^ + 8gi,292,3 - 6gi,293,4 > | = (1.4) 

J H\ J /3i 

I d^ c?/?|/9^(< 9i,2 > -6 < 9i,2 >< 92,3 > +6 < 91,293,4 >-< 9i,2 >^)| = , 

^Mi >^/3i 

(1.5) 

thus reducing to the standard ones [Pa, G, AC]. Notice the remarkable fact that the iden- 
tities of type (1.1) have recently appeared in the theory of mean field diluted ferromagnets 
[ABC]. 

In the Nishimori line [Nl, N2], 

PA'^I, (1.6) 

the identities (1.1) and (1.2) hold pointwise in P and /j, and the identities (1.4) and (1.5) 
hold when integrated over /i on arbitrary intervals. 

The plan of the paper is the following: in the next Section we define the general class 
of models for which our result apply and we set the notations. In Section 3 we state the 
theorems, which are then proved in Section 4. 
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2 Definitions 



We consider a disordered model of Ising configurations o"„ = ±1, n G A C Z'^ for some d- 
parallelepiped A of volume |A|. We denote Sa the set of all a = {cr„}neA, and |Ea| = 21^1 
In the sequel the following definitions will be used. 

1. Hamiltonian. 

For every A C Z*^ let {i7A(c)}o-eSiv be a family of 2l^l translation invariant (in 
distribution) Gaussian random variables defined, in analogy with [RU], according 
to the general representation 

Hk{(j) = - J] JxCJx (2.7) 

XcK 

where 

(Jx^\{(Ji, (2.8) 
((j0 = 0) and the J's are independent Gaussian variables with mean 

Av( Jx) = lix , (2.9) 

and variance 

M{Jx-^^xf)=^l. (2.10) 

2. Average and Covariance matrix. 
The Hamiltonian has average 

B^{a) := MH^{a)) = J] , (2.11) 



and covariance matrix 



CA(a,r) := Ky {[Hp,{a) - BK{a)][Hj,{T)) - Bk{t)]) 

= J] A^axTx . (2.12) 



XcA 

By the triangular inequality 



XcA 



for all a and by the Schwarz inequality 

|CA(a,r)| < = Y.^x (2.14) 

for all (7 and r. 

3. Thermodynamic Stability. 

The Hamiltonian (2.7) is thermodynamically stable if there exist constants c such 
that 



1 , , 
sup — 2^ |/ix| < c < oo 

sup 



A| 

^ 5] < c < oo . (2.15) 



XcA 

Thanks to the relations (2.13) and (2.14) a thermodynamically stable model fulfills 
the bound 

Ba((7) < c|A| 
CA(a,r) < c|A| (2.16) 

and has an order 1 normalized mean and covariance 

^a((7) := jxf^^^) 
CA(a,r) := ^CA(a,r). (2.17) 



4. Random partition function. 



Z{(3) := e-^^^(^) . (2.18) 



5. Random free energy. 

-(3J^{(3) := A{f3) := \nZ{/3) . (2.19) 



6. Random internal energy. 
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7. Quenched free energy. 

-(3F{(3) - A{P) := Ay {Am 



(2.21) 



8. R-product random Gibbs-Boltzmann state 



. p-/3[HA(a(l)) + -+//A(^(«')] 

"H- E H ^TmiK ■ 



o-(i),...,o-(«) 

9. Quenched equilibrium state. 

< - >:= Av(Q(-)) . (2.23) 

10. Observables. 

For any smooth bounded function G'(6a, ca) (without loss of generahty we consider 
|G| < 1 and no assumption of permutation invariance on G is made) of the mean 
and covariance matrix entries we introduce the random (with respect to < — >) 
i?-dimensional vector of elements {rrik) (called generalized magnetization) and the 
R X R matrix of elements {qk,i} (called generalized overlap) by the formula 

< G{m, q) > := Av {n{G{bA, ca))) . (2.24) 
E.g.: G(6a,ca) =6A(ai)cA(a«,a(2))cA(a(2),a(3)) 



< m,q,,2q2,3 > = Av ( J] bA{a')cA{a^'\a^'^)cA{a^'\a^'^] 

.aWM^'Kcrii) ^ ^^^^ 



(2.25) 



3 Theorems 

To state our results we introduce the random variable J' defined by 
and deform uniformly the averages iix with a parameter fi defined by 
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in such a way that 

Jx = Jx + l^l^'x ■ 
Our results can be summarized in the following theorems. 



3.1 Identities in /3-average 



For every observable G of the kind considered in the previous Section, we define 

R 



fi{P,fi) = < {mi - ruR+i) G > + 



(3.26) 



R 



R 



< ^Gqi,),-2RG ^qi^R+x + R{R + l)G QR+1,R+2 > 



k,l=l 



and 



/2(/?, ^i) = < rriR+iG > - < mi >< G > + 



(3.27) 



-P < Gqk,R+i > -{R + 1) < GqR+i,R+2 > -<G>{< qi,i > - < ?i,2 >)j 

We then have the following 

Theorem 1 The quenched equilibrium state of a thermodynamically stable Hamiltonian 
fulfills, for every observable G and every temperature interval (^2] the following iden- 
tities in the thermodynamic limit 



lim 



[ 'd/3|/i(/3,/x)| 

Jl3i 







(3.28) 



hm r\p\f2{P,pi)\ = 



(3.29) 



Remark 1 The two previous relations when applied to G{m, q) = f3qi^2 + ^1 (ind R = 2 
yields the identities mentioned in the introduction, formulae (1.1) and (1.2) - here it is 
assumed that = 1, as it happens for the Edwards Anderson model. 
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3.2 Identities in ytx)- average 
We also define 



R 



< J2 Gqi,k-2RG ^qi,R+i + R{R+l)GqR+i,R+2 > 



k,l=l 



(3.30) 



and 



^R+l 



^2(/3,Af) = -(3\^<Gqk,R+i>-{R+l)<GqR+i,R+2>-<G>{<qi,i>-<qi^2>) 

(3.31) 



,fe=i 



Theorem 2 The quenched equilibrium state of a thermodynamically stable Hamiltonian 
fulfills, for every observable G and every set [/?i,/32] x [111,112] the following identities in 
the thermodynamic limit 

lim / dn / df3\gi(i3,n)\ = (3.32) 

/•P2 

lim / diJ, / di3\g2(f3,ii)\ = (3.33) 

Remark 2 Two similar families of identities have been proved in [CGi2] for the centered 
case Av {Jx) = for all X. This theorem generalizes the old result and reduces to it when 
the observable G doesn't depend on m. 

3.3 Identities pointwise 

Theorem 3 In the space of parameters {/ix, ^x)xeA there exists a region called the 
Nishimori manifold 

f^x = PAj, (3.34) 
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where the identities of Theorem (1) hold pointwise, namely 

lim = 

A yv.d 



(3.35) 



lim /2(/3,/i) = (3.36) 

4 Proofs 

Theorem 1 is proved in the lemmas of subsections 4.1 and 4.2. The proof uses only 
elementary methods like concentration and classical inequalities. Along the same line 
Theorem 2 is proved in subsection 4.3. Theorem 3 is proved in subsection 4.4 making use 
of an exact computation on the Nishimori manifold. 

Let h{a) = |A|^^ifA(cr) denote the Hamiltonian per particle. We consider the quantity 

R 

J2 {< K(^^^^) G> - < h{a^^^) ><G>]= AiG + A2G (4.37) 
1=1 

where 

R 

AiG = {^^ HKa^^^) G] - n[h{a^^^)]n[G])} (4.38) 
1=1 

R 

A2G = J2 i^[h{a^^^)]fl[G]) - Av {fl[h{a^^^)]) Av {fl[G])} (4.39) 
1=1 

We are going to show that both AiG and A2G vanish (in P average) in the thermodynamic 
limit. This implies, by a simple application of integration by parts, the relations (3.28) 
and (3.29). 

4.1 Stochastic Stability Bound 

We follow the method of stochastic stability as developed in [CGi] . 
Lemma 4.1 For any inverse temperature interval [Pi, P2] one has 

r Av {n{h') - n^h)) d(3 < + (4.40) 
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Proof. 

The variance of the Hamiltonian per particle h{a) with respect to the Boltzmann state is 
nothing but (minus) the derivative of the average of h{a) up to a factor 1/|A|, i.e. 

Av {Q{h') - Q\h)) = -^^Av im) (4.41) 

Apphcation of integration by parts (for its general form see formula (4.50) below) yields 

= m^"^ f ^ f^xniax)] + ^ E /3A^[1 - Av {n'{ax))] 

' ' \XcA / I I XcA 

< {l + P)c, (4.42) 

where thermodynamic stability condition, Eq. (2.16), has been used in the last inequality. 
The lemma statement follows from integration of Eq. (4.41) over an arbitrary inverse 
temperature interval [^1,^2] and the use of fundamental theorem of calculus together 
with the bound (4.42). □ 

Lemma 4.2 For every bounded observable G, see definition 10 of Section 2, we have that 
for every interval [/3i,/32] in the thermodynamic limit one has 

/ lAiCl dp ^ (4.43) 

Proof. 

Prom the definition of AiG, Eq.(4.38), we have 

'■132 R 



rP2 rP2 

/ \AiG\dp < / ^ |Av(n[/i((7('))G'] -0[/i((7('))]0[G']) I d/? (4.44) 

< f'Y. ^Av ({Q[/i((7(0) G] - Q[/i((7(0)]Q[G']}') d(3 (4.45) 

< R ^J Ay {Q[h?{a)] - n'^[h{a)]) dp (4.46) 

< R^/p^^i^J^' Ay{n[h^{a)]-n^[h{a)]) dp (4.47) 

< R^M^^^^^^j^ (4.48) 
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where (4.44) follows from triangular inequality, (4.45) is obtained by applying Jensen 
inequality on the measure Av (— ), (4.46) comes from application of Schwarz inequality to 
the measure fi(— ) and boundedness of G, (4.47) is again Jensen inequality on the measure 
jp^{—)d(3 and finally (4.48) comes from lemma 4.1. □ 

Remcirk 3 The previous lemma is related to a general property of disordered systems 
which is known as stochastic stability (see [AC, CGi]). It says that the equilibrium state 
in a spin glass model is invariant under a suitable class of perturbation in all temperature 
intervals of continuity. The result presented in ([CGi]) for the case of zero average cou- 
plings holds with the absolute value outside the integral in beta with a vanishing rate of 
the inverse volume. Here instead we got the stronger result for the absolute value inside 
the integral but with a weaker vanishing rate of the square root inverse volume. 

Lemma 4.3 The following expression holds: 

R 



AiG = < {mi - niR+i) G > 



1=1 

R R 



- (3<G 



QR+1,R+2 



k,l=l 1=1 



> . (4.49) 



Proof. 

For each replica I {1 < I < R), we evaluate separately the two terms in the sum of the 
right side of Eq. (4.38) by using the integration by parts (generahzed Wick formula) for 

correlated Gaussian random variables, Xi,X2, . . . ,Xn with means Av (a;,) and covariances 
Av {{xi — Av {xi)){xj — Av {xj))), namely 

Av{Xi^{Xi,...,Xn)) = Av{xi)Av{^{xi,...,Xn)) (4.50) 

+ Av {{X, - Av {xM^j - Av (xj))) Av ( ^"^^^'g^-;' ^""^ ) ■ 

It is convenient to denote by p (R) the Gibbs-Boltzmann weight of R copies of the system 

e-^[Efc=i-ffA(aW)] 
PR ■ ■ ■ O = ^(^p , (4.51) 
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so that we have 



/5 

We obtain 



(4.52) 



\ <t(i),...,ctW / 

= Avj Yl G6A(aW)p«(<7\...,0 

+ Av( E E«^A(^'V) ^''"',-;;.--''' ) (4.53) 

\ <t(i),...,<tW r / 



= <miG> -(3 



R 



Y < Gqi,k >-R< Gqi,R+i > 



k=l 



(4.54) 



where in (4.53) we made use of the integration by parts formula and (4.54) is obtained 
by (4.52). Analogously, the other term reads 

' ' \ aW r(i),...,r(«) / 

= f E E ^ ^A(a«) {a\ n,..., tr) 

\ aW r(i),...,r(«) , 

\ (t(o t(i),...,t(«) 7 yij J 



< vfiR+i G > -(3 



R+l 



Y < Gqk,R+i > -(^+1) < ^5^^+1,^+2 > 



k=l 



(4 



Inserting the (4.54) and (4.55) in Eq. (4.38) and summing over / we obtain the expression 
(4.49). □ 

4.2 Selfaveraging Bound 

The selfaveraging of the free energy is a well established property of spin glass models. 
The vanishing of the fluctuations with respect to the disorder of the free energy can 
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be obtained either by martingales arguments [PS, CGi2] or by concentration of measure 
[T, GT2]. Here we follow the second approach. Our formulation applies to both mean field 
and finite dimensional models and, for instance, includes the non summable interactions 
in finite dimensions [KS] and the p-spin mean field model as well as the REM and GREM 
models. 

Lemma 4.4 The disorder fluctuation of the free energy satisfies the following inequality: 
for all X > 

P (1^- Av(^) I > < 2exp(^-^|p^^ (4.56) 
The free energy is then a self averaging quantity, i. e. 

V{A) = Av{A^)-Av{Af < 4:Cp^\A\ (4.57) 

Proof. Consider an s > 0. By Markov inequality, one has 

F{A-Av{A)>x} = P {exp[s(^- Av(^))] > exp(sx)} 

< Av {exp[s{A - Av {A})]) exp(-sa;) (4.58) 

To bound the generating function 

Av(exp[s(^- Av(^))]) (4.59) 

one introduces, for a parameter t e [0, 1], the following interpolating function: 

(pit) = hiAvi{exp{s Av2{l-aZ{t)})} , (4.60) 

where } and >li'2{— } denote expectation with respect to two independent copies 

Xi{a) and X2{a) of the random variable X{a), which is a centered Gaussian process with 
the same covariance as the Hamiltonian H\{a), and the partition function Z{t) is 

Indeed, since H^{a) — X{a) +Ba((7), it is immediate to verify (see definition (2.19)) that 

(/)(0) = s Av {A) , (4.62) 
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and 

0(1) = lnAv(e^-^) . (4.63) 

This implies that 

Av (exp[s(^ - Av (A))]) = e'^^^)-'^^") = e^o <t''m (4 54) 

On the other hand, the derivative with respect to t can be easily bounded. Defining 

K(t) = exp(s Av2{ln Z(t)}) (4.65) 

and 

„ -fSViXi {a)-f3VT^X2 (a) -/3Ba (a) 



one has 



Av,{K{t)sAv,{^^p\a)[^^XM-^^X,{a)\}} 
V if) — ^ ^ ~A — ( x-i (4.67) 



Av^{K{t)) 

Applying the integration by parts formula (4.50), a simple computation gives 

a ^ ^, V J J 

+ ^^;i|i^(i)^T;2|x^CA(a,a)p*(a)|| 

- Av, |x(t) Av2 |5^CA(<7,T)p*(a)p*(r))|| 



and 



Av,\^K{t) Av2\^p\a)^^^X,{a)y^ = Av,\^K{t) Av2^e{a,a)Kp\a) 

- Av, I K{t) Av2 I r)^p\a)p\T)) 



Taking the difference of the previous two expressions one finds 

... Y..,r W) CA(a, r) Av, {p\t)] Av, {p\a)]] 
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Using the thermodynamic stabihty condition (2.16), this yields 



S2 



nt)\<-\A\c (4.69) 



from which it follows 

Av (exp[s(^ - Av (A))]) < exp (4-70) 
Inserting this bound into the inequality (4.58) and optimizing over s one finally obtains 



¥ {A-Av{A)>x) < exp(^-^p^j (4.71) 

The proof of inequality (4.56) is completed by observing that one can repeat a similar com- 
putation for P (^ — Av {A) < —x). The result for the variance (4.57) is then immediately 
proved using the identity 

POD 

Av {{A - Av {A)f) =2 xF{\A- Av{A)\>x) dx (4.72) 

Jo 

□ 

Lemma 4.5 The internal energy is self averaging almost everywhere in P, i.e. defining 
u = U/\A\ and V{u) = Av (m^) — Av (m)^ it holds in the thermodynamic limit 

V{u)d(3 (4.73) 



Proof. 

The result is obtained in two steps which use general theorems of measure theory. First 
from lemma 4.4 we obtain the convergence to zero almost everywhere (in f3) of the variance 
of the internal energy, then thanks to a bound on the variance of the internal energy we 
apply the Lebesgue dominated convergence theorem which gives the lemma statement. 
The sequence of convex functions A{P)/\A\ converges a.e. (in J) to the limiting value 
a(/3) of its average and the convergence is self averaging in the sense of lemma 4.4. By 
general convexity arguments [RU] it follows that the sequence of the derivatives A'{P)/\A\ 
converges to u{(3) — a' {(3) almost everywhere in (3 and also that the convergence is self 
averaging. In fact the vanishing of the variance of a sequence of convex functions is 
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inherited, in all points in which the derivative exists (which is almost everywhere for a 
convex function), to the sequence of its derivatives (see [S, OTW]). From lemma 4.4 we 
have then 

V(u) ^ (3 - a.e. (4.74) 

In order to obtain the convergence in /3-average we use the Lebesgue dominated con- 
vergence theorem. In fact we prove that the sequence of variances of u is uniformly 
bounded (in every interval [/3i,/32]) by an integrable function of (3. A lengthy but simple 
computation which uses again integration by parts gives 

Av (W^) = Av ( JxJy^M^M ) (4.75) 
\x,ycA / 

= J2 /ixAfYJ^((Jx)^^((Jy) + XI A^Q2((7x) 

X,YCA XcA 

+ 2p /ixA^ [n((7x(7y)0((7y)+n((7x) -2n(crx)J^'(cry)] 

X,YcA 

+ (3^ A|A^Av[l-02(ax)-i^'(ay) + 6f]2(ax)0'(ay)+ 

X,YCA 

-6n{ax)^{(TY)^{(Jx(JY) + n^iaxay)] (4.76) 

from which 

V{u) < \A\-^Ay{U^) < 0^(2 + 4/3 + 14/32). (4.77) 
Prom this follows (4.73). □ 

Lemma 4.6 For every bounded observable G, see definition 10 of Section 2, we have that 
for every interval [/3i,/32] in the thermodynamic limit 



hi 



\/^2G\d(3 = (4.78) 



Proof. 

Prom the definition of A2G, Eq.(4.39), we have 



/ 



I AaGl df3 < J2 {^[h(a^^^)]n[G]) - Av {n[h(a^^^)]) Av {Q[G]) \ 
hi 1=1 

< [ \/Av (l][/i2(a«)]) - (Av(l][/i(a«)]))' d/3 (4.80) 
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< R VV{u) dp (4.81) 



Pi 



< R^(52-(5^^ j\{u) d(5 (4.82) 

where (4.79) follows from triangular inequality, (4.80) is obtained by applying Schwarz 
inequality to the measure Av (— ) and boundedness of G, (4.82) is Jensen inequality on 
the measure j^^{—)d(5. The statement (4.78) follows then using the result of the 
previous lemma. □ 

Lemma 4.7 The following expression holds: 

R 



A2G = ruR+iG > - < mi >< G >) + (4.83) 

1=1 

'R+l 

-(3R 

.k=l 



< Gqk,R+i > -{R+l) < GqR+i,R+2 > - <G>{< > - < gi,2 >) 



Proof. In order to obtain A2G we are left with the explicit evaluation of the other term 
in (4.39) which simply gives 

Av (Jl(/i(<7«))) Av (Q(G')) = ^ Av ( ^ i/A(<j«) pi {aU <G> 

= Av j ^ 6;, 



a('^)pi(a) \ <G> 



<G> 



= < mi >< G > -P < G > [< > - < gi,2 >] (4.84) 
Inserting the (4.55) and (4.84) in Eq. (4.39) and summing over I we obtain the (4.83). □ 



4.3 Vanishing fluctuations of the generahzed magnetization 
Lemma 4.8 For every interval [jii,ii2\, in the thermodynamic limit 

j dii{< > - < m >2) = (4.85) 
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Proof. The proof that the generahzed magnetization has vanishing fluctuation follows 
the strategy that has been pursued so far to control fluctuation of the internal energy. 
We have 

< > - < m >2 = Av {n{bA{af) - (0(6a((t))') (4.86) 
+ Av ((Q(6a(<7))^) - (Av mbA{a))f (4.87) 

and we observe that generalized magnetization is related to the pressure by 

The fluctuations w.r.t. the Gibbs state (r.h.s of Eq.(4.86)) are easily controlled by a 
stochastic stability argument: 

J^' dMv (nib^iaY) - mbMY) = ^ di, ^ Av (J-^n{bA{a))^ , (4.89) 

where the right hand side can be bounded, integrating by parts in // by /3~^3(//2 — A*i)- 

The fluctuations w.r.t. the disorder (Eq.(4.87)) are bounded by the same argument of 
lemma 4.5. Indeed from self-averaging of the pressure per particle 

and convexity of finite volume pressure 

one deduces that also the sequence of derivatives (4.88) is self- averaging in /x-average. 
Hence 

J'^ d^i [Av {{nib^ia))') - (Av (Q(6a(<7)))^] = d^lV (J- (^^) ) - (4.90) 
Combining together (4.89) and (4.90) complete the proof of the lemma. □ 

4.4 Nishimori manifold 

Our strategy here is to prove the vanishing of fiuctuations of the Hamiltonian per particle 
h{a). This implies, following the same reasoning of the previous sections, the vanishing 
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of correlations with a generic bounded function G using the Schwarz inequality: 



\<hG>-<h><G>\<-J<h'^>~<h >V< G'^> - <G>^ . (4.91) 

Lemma 4.9 On the Nishimori manifold 

fxx = PAl (4.92) 

the random internal energy per particle is self averaging. More precisely the following result 
holds: 

<h'^ > - <h>^ < (4.93) 

Proof. This is proved by an exphcit computation of both terms in formula (4.93). Let 
us first consider 

Av{n{Jxax)) . 
We first write the definition exphcitly: 



/OO / 1 



2 Ay // ^^e^^zJzaz 

(4.94) 

We apply the gauge transformation 

for all i e A and X C A, where Tj is a 'gauge' variable fixed to 1 or -1 at each i e X 
and Tx = Yliex This change of variables leaves the integral and sums in the above 
equation invariant. Then, only the Jy in the exponent for the Gaussian weight changes: 

Ay{Q{Jxax)) (4.95) 

(Jyry - HvfW JxCTxef^^^-^^"^ 



/OO / 

r TT 1 exp (-'A^)) exp (t'-^] ■ ^''^"^^T'" 
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Since this expression holds for any assignment of ±1 to Tj, we may sum it up over all 
possible {Ti}i and divide the result by 2l^l, 



Ay{Q{Jx(Tx)) (4.96) 



(7 

The sum over r and the sum over a in the denominator cancel each other for NL (/3 = 
/xy/Ay), and we have a simplified expression 

Av(fi(Jx(7x)) (4.97) 

- iE/In (''^^vik'^'' i-'W)) 

For given let us change the integral variable as Jy — >■ JyO'v- Then a disappears 

completely and the integral is just for the average of Jx 

Av {n{Jxax)) 

= iE/In (-^^- (-4±|i)) ..e ..-.Mi 

= Av( Jx) (4.98) 

Prom the previous computation we obtain the final result for the quenched internal 
energy on the NL: 

Ay {n{H^{a))) = Ay (q (^Jxax^j =Ay(^Jx^ =Y.^^ (^.99) 
The other term in the variance (4.93) is evaluated similarly: 
Ay {n{H^{af)) = Ay {^[^JxJy 
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= Ay ^JxJy^ 

x^y X 
= J^/Ux/^y + ^A^ (4.100) 



x,y X 

Therefore, using (4.99) and (4.100) in the expression for the variance of h one finds 



□ 
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